Draft version February 2, 2008 

Preprint typeset using LAT^^X style cmulateapj v. 11/12/01 



O 
O 

(N 

D 



PROBING DARK ENERGY USING ITS DENSITY 
INSTEAD OF ITS EQUATION OF STATE 

YuN Wang\ and Katherine Freese^ 

Draft version February 2, 2008 

ABSTRACT 

The variation of dark energy density with redshift, px{z), provides a critical clue to the nature of dark 
energy. Since px{z) depends on the dark energy equation of state wxiz) through an integral, px{z) 
can be constrained more tightly than wx{z) given the same observational data. We demonstrate this 
explicitly using current type la supernova (SN la) data [the Tonry/Barris sample], together with the 
Cosmic Microwave Background (CMB) shift parameter from CMB data (WMAP, CBI, and ACBAR), 
and the large scale structure (LSS) growth factor from 2dF galaxy survey data. We assume a flat 
universe, and use Markov Chain Monte Carlo (MCMC) technique in our analysis. We find that, while 
wx{z) extracted from current data is consistent with a cosmological constant at 68% C.L., px{z) (which 
has far smaller uncertainties) is not. Our results clearly show the advantage of using px{z), instead of 
Wxiz), to probe dark energy. 

Subject headings: cosmology: observations - distance scale - supernovae:general 
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1. INTRODUCTION 

Recent observations of type la Supernovae (Riess et al. 
1998; Perlmutter et al. 1999) indicate that the universe is 
accelerating. A fundamental quest in physics and cosmol- 
ogy is to identify the nature of the "dark energy" driving 
this acceleration. Possibilities include: (1) a cosmologi- 
cal constant, (2) a time dependent vacuum energy, or a 
scalar field known as "quintessence" that evolves dynami- 
cally with time (Freese et al. 1987; Peebles & Ratra 1988; 
Wetterich 1988; Frieman et al. 1995; Caldwell, Dave, & 
Steinhardt 1998; Zlatev, Wang, & Steinhardt 1999) ^ or (3) 
modified Friedmann equation, e.g. the Cardassian models 
(Freese & Lewis 2002; Freese 2003; Gondolo & Freese 2003; 
Wang et al. 2003), that could result as a consequence of 
our observable universe living as a 3-dimensional brane in 
a higher dimensional universe. Other proposed modifica- 
tions to the Friedmann equation include Parker & Raval 
(1999); Deffayet (2001); Bilic, Tupper, & VioUier (2002); 
Ahmed et al. (2002); Capozziello (Carloni); Carroll et al. 
(2003); Meng & Wang (2003); Puetzfeld & Chen (2004). 
The various dark energy models produce dark energy den- 
sities px{z) with different redshift dependences. Hence, 
in order to differentiate between dark energy models, it is 
important that we allow the dark energy density to be an 
arbitrary function of redshift z (Wang & Garnavich 2001; 
Wang & Lovelace 2001; Wang et al. 2003). 

A powerful probe of dark energy is type la supernovae 
(SNe la), which can be used as cosmological standard can- 
dles to measure how distance depends on redshift in our 
universe. The luminosity distance dL{z) = (1 + z)r{z), 
with the comoving distance r[z) given by 



r{z) = ci?o 



dz' 



with 

E{z) = + zf + nk{i + zf + nxpx{z)ip_ 



(1) 

I 1/2 

(2) 



where fi^ = 1 — ri„i — ^x, and px{z) is the dark energy 
density. 

The dark energy equation of state, wxiz), is related to 
px{z) as follows (Wang & Garnavich 2001): 



wx{z) = ^{l + z) -- 

3 Px{z) 



1, 



so that 



Pxjz) 
px{0) 



exp 



3[l + wxiz)] 
1 + z 



(3) 



(4) 



One can see that it is easier to extract px{z) from the data 
than to extract wx{z). To obtain the dark energy density 
directly, one need only take a single derivative of the lu- 
minosity distance, whereas to extract wx(z), one needs 
to take a second derivative as well; from Eq.(3) one can 
see that wxiz) is on the same footing as p'x{z). Specifi- 
cally, Wang & Garnavich (2001) argued that px{z) should 
be preferred since it suffers less from the smearing effect 
(due to the multiple integrals that relate wx{z) to dL{z)) 
that makes constraining wx{z) extremely difficult (Maor, 
Brustein, & Steinhardt 2001; Barger & Marfatia 2001). 
Tegmark (2002) came to the same conclusion. However, 
researchers have generally chosen to parametrize dark en- 
ergy using its equation of state wx{z)- Some have used 
H{z) = HqE{z) (for example, see Kujat et al. (2002); Daly 
& Djorgovski (2003); Nesseris & Perivolaropoulos (2004), 
and references therein), which is similar to px{z), but mea- 
surements of which are not as straightforward to interpret, 
since E{z) depends on (see Eq.(2)). 

In this paper, we explicitly demonstrate the advan- 
tage of using px{z), instead of wx{z), to probe dark en- 
ergy. Sec.2 contains a comparison of wx{z) and px{z) 
parametrizations using current SN la, CMB, and LSS 
data. We give a recipe for parametrizing dark energy using 
px{z) in Sec. 3. Sec. 4 contains a summary and discussions. 
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and dashed lines indicate the 68.3% and 95% confidence 
levels respectively (Wang & Mukherjee 2004). The other 
simultaneously estimated parameters (mean, 68.3% and 
95% confidence ranges) are: fi„ = .33[.27, .39] [.22, .46], 
h = .660[.644, .673][.630, .688]. One can see that the un- 
certainties in Fig. lb on px{z) obtained from the data are 
smaller than those on Wx{z) obtained from the data. We 
see that the time dependence of the dark energy density 
deviates from a constant at 68.3% C.L. (a similar state- 
ment could not be made from the wx{z) reconstruction). 
With more data in the future, the statistical significance 
of this discrepancy will become more clear. 

For comparison, in Fig. 2 we have also plotted px{z) ob- 
tained in a more indirect way: by first obtaining wx{z) 
from the data, as described above, and then integrating 
over redshift as in Eq.(4). Clearly the uncertainties ob- 
tained in this way are far larger than if one obtains the 
dark energy directly from the data. While the results from 
wx{z) parametrization and px{z) parametrization (Wang 
& Mukherjee 2004) are consistent with one another, the 
wx {z) parametrization results have uncertainties that are 
several times larger. One can also obtain wxiz) indirectly 
from px{z) in Fig. lb (similar to what was done by Alam 
et al. (2003)). However, doing so would require taking the 
derivative of the polynomial used in the interpolation, thus 
making the resultant Wx{z) dependent on the interpola- 
tion technique used. 

Our main result is that one can learn more information 
by reconstructing px{z) rather than wx{z) from the data. 
At 95% C.L., both the wx{z) and px{z) reconstructions 
are consistent with a cosmological constant. However at 
68.3% C.L., the px{z) reconstruction has smaller imccr- 
tainties and hence shows more information than the wx{z) 
reconstruction: the px {z) reconstruction is not consistent 
with a time-independent dark energy. Even with the px {z) 
parametrization, a significant number of SNe la at z > 1 
from a deep SN survey on a dedicated telescope (Wang 
2000a) will be required to place robust constraints on the 
time-dependence of px{z). 

3. A RECIPE FOR PARAMETRIZING DARK ENERGY 
USING ITS DENSITY 

For the convenient application of our methodology by 

others, we now present a recipe for parametrizing dark 
energy using px{z) as an arbitrary continuous function. 

(1) Choose the number of redshift bins, N (the number 
of parameters for px{z)). N needs to be sufficiently large 
to probe the time- variation of px{z). However, if N is too 
large, the uncertainties on all the estimated parameters 
will increase, leading to less stringent constraints. N = 2 
is appropriate for current (sparse) data. 

(2) The values of the dimensionless dark energy den- 
sity fi = p{zi)/px{0) {i = 1,2,...,N) are the independent 
variables to be estimated from data. Note that zjsi = Zmax 
(the maximum redshift of SNe la in the data). 

* Here is only used to move around efficiently in the entire parameter space (based on entropy considerations), such that for sufficent 
sampling, the resultant parameter distributions converge to the true probability distribution functions (pdf's). This leads to smooth pdf's since 
they receive contributions from all MCMC samples. 

'' This is the largest number of values cue can get out of the current (sparse) data. 

^ Flux averaging has been performed to reduce the bias in estimated parameters due to weak gravitational lensing. See Wang (2000b) and 
Wang & Mukherjee (2004) for details. 

The uncertainty on h axe statistical error only, not including the contribution from the much larger SN la absolute magnitude error of 
a'^* ~ 0.05. (Wang & Mukherjee 2004) 



2. DARK ENERGY EQUATION OF STATE VERSUS 
DARK ENERGY DENSITY 

;^From SN, CMB, and LSS data, we independently re- 
construct first the dark energy equation of state and then 
the dark energy density directly. We use (1) current 
SN la data from Tonry et al. 2003 and Harris et al. 
2003 [the Tonry/Barris sample], (2) the Cosmic Microwave 
Backgroimd (CMB) shift parameter (Bond, Efstathiou, 
& Tegmark 1997) from CMB data [WMAP (Bennett et 
al. 2003; Spergel et al. 2003), CBI (Pearson et al. 2003), 
and ACBAR (Kuo et al. 2002)], and (3) the large scale 
structure (LSS) growth factor from 2dF (Percival et al. 
2002; Verde et al. 2002; Hawkins et al. 2003) galaxy sur- 
vey data. (Wang & Mukherjee 2004) 

We parametrize the data first in terms of the dark en- 
ergy equation of state, and then in terms of the dark en- 
ergy density, to see which parametrization produces the 
reconstruction with the least uncertainty. 

To find the function wx (z) which best fits the data, we 
consider a five-dimensional parameter space: the function 
evaluated at three discrete redshift intervals, the value of 
i^m, and the value of the Hubble constant. We assume a 
flat universe. It is our goal to find a set of values for these 
five parameters that fits the data. We use the Markov 
Chain Monte Carlo (MCMC) technique (Neal 1993; Lewis 
& Bridle 2002), which selects randomly from the 5D pa- 
rameter space, evaluates x^, and create a large number of 
sets of parameter values (each set is a MCMC sample); 
we use 10® MCMC samples ^. For the current data, we 
consider the function wxiz) at three redshift values: z=0, 
Zmax/ '2', and Zmax (where Zmax is the maxmimum redshift 
of SNe la)^. We find the values (that fit the data) at these 
points and interpolate at all intermediate redshifts. The 
parameters estimated from data are wx{0), wxizmaxf^), 
wxizmax), ^m, and & dimeusiouless Hubble constant h. 

Fig. la shows the Wxiz) reconstructed from 192 SNe 
la from the Tonry/Barris sample, ® combined with CMB 
(shift parameter TZq = 1.716±0.062) and LSS data (growth 
parameter /o = f(z = 0.15) = 0.51 ± 0.11)(Wang & 
Mukherjee 2004). The regions inside the solid and dashed 
lines correspond to 68.3% and 95% confidence levels re- 
spectively; the 68.3% confidence level (C.L.) region is also 
shaded. The circles indicate the mean values at the three 
redshift points, z=0, Zmaxf^^ and z„iax- The other simul- 
taneously estimated parameters (mean, 68.3% and 95% 
confidence ranges) are: flm = ■39[.29, .50][.21, .57] and 
h = .658[.642, .674][.627, .689].'^ Clearly, the equation of 
state is consistent with a constant wx{z) = — 1 for all red- 
shifts at 95% confidence level (C.L.). At 68.3% C.L., it is 
consistent with a constant for < z < 0.5 and marginally 
consistent with wx{z) = —1 for 0.5 < z < 1. 

Fig. lb shows the pxiz) directly reconstructed from the 
same data as Fig. 1. The same technique of discretiz- 
ing the function px{z) has been used. The solid lines 
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(3) Parametrize px{z)/ px{^) as a continuous function, 
given by interpolating its amplitudes at equally spaced 
z values in the redshift range covered by SN la data 
(0 < z < Zraax), and a constant at larger z {z > z^axi 
where px^z) is only weakly constrained by CMB data). 
The results should not be sensitive to the interpolation 
method used. Polynomial interpolation was used in Wang 
& Mukherjee (2004). For A/" = 2, this gives 



Px{z) 
Px(0) 



= 1 + (4/i - /2 - 3) 



z 2z^ 



(5) 

where /i = p{zmax/'^)/ Px{Q)-, and = p{zmax)/px{0)- 

(4) Use Eq.(5) or its equivalent (if the interpolation 
method or N differs) in all equations where the factor E{z) 
from Eq.(2) appears. 

Caution: It is important to note that we are using a 
polynomial to interpolate px {z) between equally spaced z 
values; the independent variables are the values of px{z) 
at these z values, as in Eq.(5). In this case, the errors 
on the reconstructed px{z) are tied to how the quality of 
the data varies with z (sparse data lead to large errors). 
Changing the interpolation method from polynomial in- 
terpolation to a different method should have negligible 
effect on the reconstructed px{z). On the other hand, if 
a polynomial is used as a global fit function with its co- 
efficients being the independent variables; the errors on 
the reconstructed px{z) will not correlate with how the 
quality of data varies with z. 

4. SUMMARY AND DISCUSSION 

The critical first step in solving the mystery of dark en- 
ergy is to determine whether the dark energy density px {z) 
varies with time. (Wang & Garnavich 2001) A definitive an- 



swer to this question can have profound implications for 
particle physics and cosmology. 

Our main result is that one can learn more information 
by reconstructing px{z) rather than Wx{z) from the data. 
The two quantities are related by an integral, which in 
the case of wx{z) smears out much of the information one 
could otherwise learn. We show this explicitly by using a 
combination of SN la data from the Tonry/Barris sample 
as well as CMB (WMAP, CBI, and ACBAR) and large 
scale structure (2dF) data. At 95% CL, both the wx{z) 
and px{z) reconstructions are consistent with a cosmo- 
logical constant. However at 68% CL, the px{z) recon- 
struction has smaller uncertainties and hence shows infor- 
mation that the wxiz) reconstruction cannot: the px{z) 
reconstruction is not consistent with a time-independent 
dark energy, and the dark energy density appears to be 
increasing with redshift. Future data will be required to 
resolve this question. 

We have shown definitively the advantage of the px{z) 
parametrization over the Wx{z) parametrization in deter- 
mining the time- variation of px{z). To help others apply 
the px{z) parametrization, we have given a recipe for us- 
ing the px{z) parametrization in data analysis to probe 
dark energy (see Sec. 3). Our methodology should be very 
useful in all data analysis aiming at unraveling the nature 
of dark energy. 
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Fig. 1. — (a) The wx{z) reconstructed from 192 SNe la from the Tonry/Barris sample, flux-averaged, and combined with CMB and LSS 
data. The regions inside the solid and dashed lines correspond to 68.3% and 95% confidence levels respectively; the 68.3% confidence level 
region is also shaded. Circles indicate the mean values of the regions, (b) The px(^) reconstructed from the same data as Fig. 1(a), with the 
same shading and line types (Wang & Mukherjee 2004). Whereas the wx(^z) reconstruction is consistent at 68.3% C.L. with a cosmological 
constant, the px{z) reconstruction is not. 
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Fig. 2. — The px{^) reconstructed by taking the integral in Eq.(4) of the wx{z) plotted in Fig. 1(a). Again, we use the same data, with 
the same shading and line types. Clearly the uncertainties in this method are much greater than if one obtains px (z) directly from the data 
as in Fig. 1(b). 



